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COMPARISON  WITH  OBLIQUE  SOUNDER  DATA  OF  HICH - LATITUDE  HF  PROPAGATION  — 1 - 

PREDICTIONS  FROM  "RADAR  C*  AND  "AMBCOM"  COMPUTER  PROCRAMS 

Mat 

Hlkhll  Dave* 

Oca an  and  Ataoapheric  Sclancaa  Dlvlaton  (l 

Naval  Ocaan  Systems  Cancar  1. 

San  Diego ,  CA  92152-5000  fl. 

ABSTRACT  *  * 

A  atudy  la  dona  using  two  HF  propagation  prediction  program*  •  "RADAR  C* 
and  "AMBCOM"  •  Co  dateralna  how  wall  ehay  predict  aadlan  value*  of  oblique 
aoundar  data  of  aaxlaua  obaarved  fraquanelaa  (MOF)  aC  high  latitude*.  The 
Bain  difference*  between  RADAR  C  and  AMBCOM  are  the  Inclualon  In  the  latter  of 
hlgh-latltuda  lonoaphere  and  auroral  abaorptlon  aodela,  a a  wall  aa  a  more 
aophlatlcacad  and  accurate  ray-tracing  acheae.  The  data  used  for  coaparlaon 
are  taken  froa  Reference  [1]  for  the  Vinnlpeg-Reaolute  Bay  path  In  the  year 
1959  (alao  discussed  by  Petrie  and  Uarren  [2])  and  froa  Folkestad  [3]  for  the 
Andoys-Ft.  Honaouth  and  Andoya -Col lege  paths  In  1964.  The  data  for  the 
Vlnnlpag-Resolute  Bay  corresponds  to  high  sunspot  nuaber,  while  the  others 
correspond  to  low  sunspot  nuaber.  Hence,  this  study  provides  lnforaatlon  on 
the  perforaance  of  the  two  programs  for  various  hlgh-latltude  paths  at  both 
high  and  low  sunspot  nuaber. 

AMBCOM  was  found  to  give  generally  better  agreeaent  with  Che  above  data 
than  did  RADAR  C.  Coaparlaon  of  details  of  nodal  predictions  froa  the  two  coa- 
puter  programs  for  the  above  data-base  la  used  to  fora  an  understanding  of 
this  improvement  In  prediction  capability. 

INTRODUCTION 

This  paper  begins  with  a  summary  of  the  differences  between  the  basic 
Ionospheric  aodals  and  raytraclng  assumptions  Bade  in  constructing  the  RADAR  C 
and  AMBCOM  prograaa.  User  options  selected  for  this  study  are  discussed  in 
Section  2.  In  Section  3,  coaparlaon  of  predictions  froa  the  two  programs  with 
available  oblique  sounder  data  Is  presented  with  appropriate  explanation.  In 
Section  4,  the  coaparlsons  with  data  are  discussed  In  terms  of  what  they 
reveal  about  the  significances  of  the  differences  becveen  the  two  programs, 
and  conclusions  are  formed  regarding  the  apparent  reasons  for  Improved  predic¬ 
tive  capability  of  AMBCOM  over  RADAR  C.  In  Che  final  section,  suggestions  for 
directions  In  future  work  cowards  improving  HF  propagation  prediction  In  hlgh- 
latltude  regions  are  made. 

SECTION  1  -  SOME  BASIC  ASSUMPTIONS  USED  IN  THE  DEVELOPMENT  OF  RADAR  C 
AND  AMBCOM  COMPUTER  PROCRAMS 

The  developmental  histories  of  RADAR  C  and  AMBCOM  are  different,  and  this 
fact  accounts  for  soae  of  the  differences  becveen  Che  two  program* .  RADAR  C 
was  developed  Co  predict  perforaance  of  over-the-horizon  radars  (Headrick,  at. 
al.  [4],  Lucas,  at.  al.  (5]).  Thus,  RADAR  C  has  only  a  coverage  option,  not  a 
point- co-polnc  or  "hoaing”  option.  The  propagation  nodal  Is  based  on  virtual 
geometry  and  Is  essentially  ch*  saae  as  chat  of  ITS  -78  (Barghausen,  at.  al. 
(6) )  and  IONCAP  (Tatars,  at.  al..  [7]).  AMBCOM  was  derived  froa  the  NUCOM 
prograa  developed  at  SRI  International  (under  the  sponsorship  of  Che  Defense 
Atomic  Support  Agency,  DASA,  and  its  successor,  Che  Defense  Nuclear  Agency, 
DNA) .  The  purpose  of  NUCOM  is  to  predlcc  the  effects  of  a  nuclear  disturbance 
on  Ionospheric  communication  channels  (Nielson,  et.  al.  [8]),  and  as  a  part  of 
thla'_ohjecclv* ,  AMBCOM  was  developed- to  predict  HF  propagation  In  an  undis¬ 
turbed,  or  ambient.  Ionosphere.  The  raytraclng  scheme  In  NUCOM/AMBCOM  was 
developed  specifically  to  permit  Che  treatment  of  a  non-hor lzontally 
stratified  ionosphere  In  the  direction  of  propagation  (1.*.,  It  Includes 
modelling  of  longitudinal,  but  not  transverse,  tilts,  so  that  propagation  Is 
along  the  great  circle  path).  AMBCOM  has  boch  coverage  and  point-to-point  op¬ 
tions  . 
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(1)  Determination  of  Ionospheric  puiMtiti  lo  dono  using “the  io-c>ll«d  ESSA 
•blue  dock*  coefficients  In  AMBCOH  with  high- loti  tudo  nod  If  lcatlons  Introduced 
by  SRI  International  (Hatfield  (10]).  RADAR  C  uses  the  unmodified  ESSA  ‘red 
deck*  coefficients. 

(2)  AMBCOH  chooses  up  to  41  control  points  (depending  on  path  length)  to 
determine  local  Ionospheric  paraaaters  such  as  critical  frequencies,  whereas 
RADAR  C  has  a  maximum  of  4  control  points  available  for  the  user  to  Input. 

(3)  AHBCOH  node  la  the  Ionosphere  with  three  parabolic  layers  of  electron  den¬ 
sity  as  a  function  of  height,  and  uses  a  seal - analytic ,  two -dimensional 
ray c racing  scheme  based  on  a  method  due  to  Rift  and  Fooks  (Nielson  (11]).  The 
physical  bases  of  this  scheme  are  the  geometric  optics  solution  to  the  wave 
equation  and  Fermat’s  principle  of  minimum  phase  (Kelso  [12]).  RADAR  C  uses 
vertical  lonograms  computed  froa  a  similar  Ionospheric  model  as  AMBCOH,  and 
converts  to  oblique  propagation  using  Martyn’s  theorem  (Davies,  (13]).  The 
scheme  used  in  AMBCOH  permits  consideration  of  continuous  Ionospheric 
gradients  along  tha  direction  of  propagation,  whereas  tha  RADAR  C  scheme  as¬ 
sumes  horizontal  stratification  of  the  Ionosphere  at  each  reflection  point. 

An  added  difference,  not  explicitly  noted  In  T4rt>le  1,  Is  the  fact  that 
AMBCOH  Is  bettar  capable  than  RADAR  C  of  considering  composite  modes  Involving 
reflections  from  tha  E,  FI,  and  F2  layers.  Including  topside  reflections  off 
of  the  lower  layers  (K-modas),  as  well  as  chordal  or  perigee  modes  (l.e.,  rays 
which  do  not  Intersect  the  earth  between  layer  reflections),  as  possible  modes 
of  propagation.  This  improved  capability  of  AMBCOH  is  due  to  Its  more  ac¬ 
curate  raytrace  method,  (a.  g. ,  In  AMBCOH  it  Is  not  assumed  that  Che  angle  of 
Incidence  to  a  layer  equals  the  angle  of  reflection).  Although  RADAR  C  Is 
also  capable  of  considering  composite  modes,  the  assumption  of  horizontal 
stratification  prevents  Che  consideration  of  tilts  and  chordal  modes  by  this 
program,  as  a  result  of  which  the  majority  of  nodes  found  by  RADAR  C  turn  out 
to  be  simple  modes  (all  reflections  being  off  of  Che  same  layer).  In  summary, 
the  treatment  of  modes  in  AMBCOH  Is  closer  to  physical  reallcy  chan  chat  In 
RADAR  C. 

It  should  be  noted  here  that  neither  of  Che  programs  Is  designed  ex¬ 
pressly  for  predicting  maximum  usable  frequency  (MUF)  for  a  given  model 
ionosphere,  so  Chat  the  program  output  has  to  be  interpreted  to  estimate  a 
MUF.  For  this  study,  since  RADAR  C  does  noc  have  a  point-to-point  option.  Its 
output  for  a  given  condition  Is  scanned  for  the  maximum  frequency  whose  ground 
ranges  (for  some  takeoff  angle)  bracket  the  receiver,  this  being  Interpreted 
as  the  HUF.  In  Interpreting  AMBCOH  output,  the  point-to-point  option  Is 
chosen,  and  It  Is  assumed  that  all  modes  which  reach  the  ground  within  100 
kilometers  ground  distance  of  Che  receiver,  or  all  chordal  modes  which  reach 
less  than  °0  kilometers  heleht  above  rhe  receiver  ire  deferrable  -orle*  The 
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ran^«*  In  thaaa  acceptance  ertcarla  ar*  somewhat  broader  than  normally  uaad 
(20  km  height  being  a  Bora  common  limit  for  chordal  aodaa,  for  example),  but 
ar*  ballavad  to  b*  rapraaaneaclv*  of  tha  rang*  of  dlacancoa  from  which  aodaa 
can  ba  dacactad,  conaldarlng  tha  accuracy  of  eh*  lonoapharlc  modal  and 
raycraclng  achaaa,  and  broadening  of  tha  baaa.  Tha  raaulta  of  thla  aeudy, 
Judgad  by  examining  tha  output,  ar*  not  highly  aanaltlv*  to  eh*  cholc*  of  tha 
above  number*, 

SECTION  2  •  USER  OPTIONS  TAKEN  IN  PERFORMING  THE  STUDY 

Soma  of  the  uaar  option*  available  In  the  two  program*  are  of  relevance 
to  thla  aeudy,  hone*  are  dlacuaaed  below: 

(1)  Four  control  point*  ar*  uaad  In  thla  aee  of  RADAR  C  run* ,  approximately 
uniformly  apaced  on  tha  graat  circle  path  between  tranamlttar  and  receiver. 

(2)  Sporadic  E  mode*  are  not  conaldarad  In  thla  study  (I0PES  -  0  In  AHBCOK, 
MAXMOD  -  1  In  RADAR  C) . 

(3)  12  month  running  average*  of  monehly  median  aunapot  number  are  uaad. 
Monthly  median  valuaa  of  magnetic  Index  R.  are  uaad  In  AMRCOM.  (RADAR  C  doe* 
not  uaa  Kp  ) . 

(4)  In  order  to  minimize  eh*  amount  of  computing  time  without  a  graat 
•acrlfle*  In  accuracy,  only  Integer  value*  of  frequencies  In  eh*  rang*  of  1  • 
30  MHz  are  Input  for  study  In  these  programs.  Thu*  tha  predicted  maximum 
usable  frequencies  have  up  to  0.5  MHz  systematic  bias  on  the  low  aid*,  alnc* 
tha  actual  MUF  would  be  leas  than  the  lowest  (Integer)  frequency  for  which  no 
propagation  la  predicted  by  the  rayerac*  scheme,  but  possibly  higher  than  the 
highest  one  found  supported. 

(3)  Tha  *hlgh-ray*  calculation  option  In  AKBCOM  Is  chosen  (HIRAY(l)  “  0.4), 
permitting  Identification  of  possible  hlgh-angle  rays  on  a  given  path. 

(6)  In  AMBC0M,  take-off  angles  from  0  to  45  degrees  are  considered,  with  one 
degree  Increments  between  angles.  In  RADAR  C,  two  degree  increments  in 
takeoff  angle  were  used,  as  It  was  merely  necessary  to  bracket  the  rang*  to 
tha  receive  alts. 

A  map  showing  the  paths  studied  la  given  In  Figure  1,  and  a  •unary  of 
geographical,  temporal,  and  solar  parameters  pertaining  to  the  data  la  given 
In  Table  2.  V*  note  that  this  study  Includes  a  short  path  for  which  most  of 
the  modes  should  be  1-hop,  and  two  Intermediate  length  paths  for  which  com¬ 
posite  mode  propagation  can  be  Important.  The  short  path  data  Is  at  high 
sunspot  number,  while  the  longer  path  data  Is  at  low  sunspot  number.  On*  of 
the  longer  paths  can  be  considered  a  trans -auroral  path  while  the  other  can  be 
considered  a  trans-polar  path  (Folksstad  [3]). 
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Iln  Figures  2a- 12b  art  p relented  tha  MUF  predictions  deduced  from  the 
two  programs  for  each  path-month  studied,  along  with  an  identification  (below 
the  universal  time,  or  UT,  axis)  of  the  mode  which  determines  the  MUF.  and  its 
:  corresponding  total  path  loss  in  dB.  for  every  two  hours  of  UT.  The  results  for 
the  two  programs  are  arranged  side-by-side,  the  figure  numbered  with  "a" 
j  corresponding  to  RADAR  C  predictions,  and  that  numbered  with  "b“  correspond- 
j  ing  to  AM B COM  predictions.  The  notation  for  modes  used  is  explained  In  Davies 
'  [13].  A  minus  sign  indicates  a  chordal  ray.  and  a  "V  indicates  topside  reflection. 
Thus,  ”£  -F2"  indicates  a  2-hop  chordal  (or  perigee)  ray  which  reflects  off  of  the 
.  £  layer,  intersects  the  earth,  then  reflects  off  of  the  F2  layer,  reaching  the  re¬ 
ceive  site  at  an  altitude  of  not  more  than  90  km  (c.f.  Section  1).  likewise.  *T2 
if  1  FZ"  Indicates  a  mode  which  reflects  off  of  the  F2  layer,  then  off  of  the  top¬ 
side  of  the  FI  layer,  then  again  off  the  F2  layer.  An  If  refers  to  a  high- angle 
mode.  _  i  — _ 
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(3)  On  the  Andoya-College  path  (Figures  9-12),  ANBC0M  shows  a  significant  im¬ 
provement  In  MUT  prediction  over  RADAR  C  (which  Is  generally  5-10  MHz  too 
low)  ,  with  composite  and  chordal  aodes  playing  an  laportant  part  in  determin¬ 
ing  the  KUF.  There  are  several  cases  in  which  2-  and  3-hop  aodes  involving  a 
coablnatlon  of  E  and  F2  layer  reflections,  as  well  as  chordal  aodea,  determine 
the  KUF.  Although  AHBC0M  Is  a  significant  Improvement  over  RADAR  C  In  this 
case,  there  ts  room  for  aore  improvement,  as  AKBC0H  Is  still  biased  on  the  low 
aide  of  observed  median  MUF’s. 
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(4)  As  a  by-produce  of  the  face  ehae  AKBCOM  generally  flnda  higher  MUF  value* 
(which  ara  closer  to  eh*  observed  values)  Chan  RADAR  C,  AMSCOM  also  shows 
lower  pach  losses  for  chess  higher  frequencies,  so  thsc  the  required  power  on 
Certain  paths  >ey  be  slgniflcancly  lower  Chan  Chat  predicted  by  RADAR  C. 

For  added  Insight,  a  comparison  of  the  calculated  values  of  E-  and  F2- 
layer  critical  frequencies  and  heights  of  the  lsyer  maxima  for  the  two 
programs  on  the  Andoya-Collaga  path  at  6  hour  Intervals  Is  plotted  In  Figure 
13.  (E-layer  maximum  height  1*  a  constant  130  km  In  RADAR  C  and  113  km  In 
AMBCOM)  .  This  figure  shows  ehat  Che  Ionospheric  parameter*  on  the  Andoya- 
College  path  differ  for  the  most  part  by  only  a  few  percent  beewesn  the  two 
programs  so  ehat  on*  can  conclude  ehat  the  radical  Improvement  In  MUF  predtc- 


tion  of  AHBCOM  ovor  RADAR  C  for  thl*  p«ch  It  not  duo  aalnly  to  tho  voluoo  of 
lonoaphorlc  poroaotoro  ut'd.  Rothor.  boood  on  tho  typos  of  »odos  found  to 
cnnatltuta  tha  KUT  in  AHBCOM.  It  i.  to  b.  eoncludad  thaV  th.  .or.  aeeurata 
*nd  physically  .or.  r.all.tlc  raytracing  in  AHBCOM.  eoabinad  with  unv  .or. 
control  points  than  uaad  in  RADAR  C,  ar.  th.  Min  caua.a  for  th.  significant 
l.proy.a.nt  in  HUT  pr. diction  in  AHBCOM  co.par.d  to  RADAR  C  Th... 
capablllti.a  (c.f.  th.  diacuaaion,  it*.  (3)  of  S.ction  1)  allow  for  ior.  »c- 
curata  tr.at.ant  of  tilt,  and  co.po.ita  and  chordal  aod.a,  which  ia  not 
poaaibl.  in  RADAR  C,  owing  to  th.  aaauaption  of  Martyn'a  th.oraa.  Thia  conclu- 
. ion  la  conalat.nt  with  th.  atudy  by  Raul  {14)  of  th.  laportince  Of  horirontal 
gradl.nta  in  al.ctron  danaity  in  th.  ionoaph.r.  .van  at  aid- laclc-ud.a 
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A  alailar  coaparlaon  of  lonoapharic  paraaetara  for  the  evo  program  for 
tha  Crana-auroral  path  (Andoya-Ft.  Monaouth)  la  ahovn  In  Figure  14.  Thla 
figure  aupporta  tha  conclualon  that  for  thla  path.  In  addition  to  tha  affacta 
of  laprovad  raytracing,  dlffarancoa  In  lonoapharic  aodalllng  (aapaclallr  for 

tha  E- layer  critical  frequency),  era  alao  algniflcant  cauaaa  for  tha  laprovad 
predictability  of  AMBCOM  over  RADAR  C.  _ 
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Obse  rvatlonal  ly ,  since  {•*  lonosondes  have  bten  operated  In  polar 
regions,  lc  la  Co  ba  axpaccad  that  a  program  of  varclcal  lncldanca  lonoapharlc 
critical  frequency  aaaauraaanca  In  cha  polar  region  will  Improve  our  capacity 
Co  nodal  HF  propagation  In  Chla  Important  raglon  of  Cha  world. 
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